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ABSTRACT  OF  THESIS 


A  TWO-DIMENSIONAL  BALANCED  MODEL  OF  INTERNAL  FRONTOGENESIS 
IN  GEOSTROPHIC  AND  ISENTROPIC  COORDINATES 


Here  we  seel^to  both  improve  and  simplify  the  method  by  which  upper  frontogenesis  may 
be  studied.  Using  a  two-dimensional  form  of  the  geostrop  hie  momentum  approximation 
in  geostrophic/isentropic  coordinates,  our  dynamic  model  reduces  to  a  predictive  equation 
for  the  potential  pseudo-density  (inverse  Rossby-Ertel  potential  vorticity),  with  associated 
diagnostic  equation  for  the  Bernoulli  function  from  which  the  wind  and  mass  fields  can 
be  calculated.  Ageostrophic  motions  are  implicit,  and  vertical  motions  retained  for  the 
adiabatic  case  employed,  by  this  choice  of  coordinates.  Initialization  of  the  domain  incor¬ 
porates  a  realistic  vertical  distribution  of  the  mass  field  along  with  upper/lower  boundaries 
which  are  either  isobaric/isentropic  or  constant  potential  vorticity  surfaces.  Vertical  wind 
shears  such  as  are  commonly  associated  with  baroclinic  waves  are  idealized  and  act  as  the 
forcing  mechanism  for  frontogenesis. 

Major  model  results  include  the  formation  of  upper  fronts  with  associated  wind  and 
thermal  fields  which,  when  viewed  together,  are  well  correlated  with  observations  of  these 
parameters  in  terms  of  magnitude  and  gradient  as  well  as  their  proximity  to  one  another. 
Well-defined  folding  of  the  dynamic  tropopause  occurs  in  geostrophic  space;  thus,  unlike 
previous  balanced  models,  the  transformation  back  to  physical  space  is  not  required  in 
order  to  produce  the  desired  results.  However,  we-show-that^performing  the  coordinate 
transformation  enhances  the  realism  of  the  results. 

Our  findings  document  the  applicability  of  the  geostrophic  momentum  approximation 
to  non-curved  flows  containing  high  relative  vorticity,  and  the  simplicity  of  the  dynamics 
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when  applied  to  these  particular  coordinates.  Further,  we  find  that  generalizing  the  forcing 
mechanism  can  also  produce  noteworthy  aspects  of  internal  frontogenesis;  this  was  accom¬ 
plished  by  extracting  an  essential  feature  of  baroclinic  waves — vertical  wind  shearing — and 
incorporating  it  first  into  the  classical  deformation  field  and  later  applying  it  to  the  zonal 
flow  directly.  The  results  suggest  more  can  be  learned  about  frontogenesis  by  application 
of  balanced  two-dimensional  theory  even  when  relatively  extreme  simplifications,  such  as 
adiabatic,  frictionless  flow,  are  made. 
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Chapter  1 


INTRODUCTION 

The  study  of  internal  (upper-level)  fronts  has  been  an  active  field  of  research  since 
the  late  1920s,  when  upper-air  instruments  first  became  operational.  Some  of  the  salient 
structures  of  these  fronts — features  such  as  static  stability  in  the  frontal  layer,  trans¬ 
verse  horizontal  temperature  gradients,  and  the  nearly  isothermal  layer  in  the  region  now 
referred  to  as  the  tropopause  and  lower  stratosphere — were  observed  early  on.  The  promi¬ 
nent  interpretation  which  grew  out  of  these  observations  was  that  of  the  Bergen  School — a 
continuous  frontal  theory — which  evolved  into  a  “deep”  polar  frontal  model  where  fronts 
were  thought  to  extend  from  the  surface  to  the  tropopause.  Not  surprisingly,  both  the 
avenues  and  roadblocks  to  knowledge  of  upper-level  fronts  have  paralleled  the  advance  of 
the  instrumentation  necessary  to  measure  their  characteristics;  for  example,  the  Bergen 
School  theory  fell  into  disuse,  at  least  initially,  primarily  due  to  the  lack  of  supporting 
observational  evidence. 

One  important  advance,  advocated  by  Reed  and  Danielson  (1959),  was  that  fronts 
were  no  longer  viewed  as  regions  separating  polar  from  tropical  tropospheric  air.  In  this 
model,  upper  fronts  are  perceived  as  a  secondary  phenomenon  of  a  tropopause  folding  event 
through  tilting  effects,  and  that  these  fronts  may  be  defined  as  separating  tropospheric 
and  stratospheric  air  instead  of  separating  two  distinct  tropospheric  air  masses.  Two  im¬ 
plications  of  this  model  are  that  we  can  specifically  account  for  stratospheric-tropospheric 
mass  exchange,  a  result  which  was  later  confirmed  (e.g.  Briggs  et  at. ,  1963),  and  that 
upper  fronts  may  form  independently  of  surface  fronts.  Observational  studies  by  Shapiro 
et  al.  (1987)  indicate  the  upper  low  pressure  vortex  and  jet  stream  associated  with  a 
folding  tropopause  occur  in  arctic  weather  systems  as  well  as  the  more  commonly  referred 


to  mid-latitude  and  subtropical  systems.  This  data  led  to  their  proposal  of  a  new  model 
for  the  meridional  structure  of  the  tropopause  and  primary  wind  systems — referred  to  as  a 
three-fold  model — wherein  there  exists  stratospheric-tropospheric  exchange  through  these 
folds  accompanied  by  subtropical,  polar,  and  arctic  jetstreams,  respectively. 

Tropopause  folding  as  a  mechanism  for  upper  frontogenesis  is  not  inconsistent  with 
the  modern  view,  which  considers  such  frontogenesis  to  be  related  to  the  development 
of  baroclinic  waves.  Kinematically,  however,  baroclinic  waves  with  their  associated  wind 
shear  patterns  result  in  a  complex  process  wherein  differentia]  advection  of  both  the  ther¬ 
mal  and  wind  fields  tends  to  hinder  attempts  at  understanding  any  single  aspect  of  the 
problem.  One  way  to  gain  information  about  these  processes  is  to  make  idealizations  in 
such  a  way  as  to  reproduce  frontal  characteristics  which  are  consistent  with  those  ob¬ 
served  as  growing  out  of  the  baroclinic  wave.  For  example,  we  can  isolate  frontogenesis 
by  idealizing  the  wind  and  thermal  fields  in  ways  known  to  lead  to  frontogenesis;  such 
an  approach  is  used  in  the  work  of  Hoskins  (1971)  and  Hoskins  and  Bretherton  (1972) 
where  a  geostrophic  deformation  field  acts  on  a  weak  initial  temperature  gradient.  This 
procedure  produces  realistic  scale  contraction  wherein  the  fluid  concentrates  into  narrow 
zones,  resulting  in  cross-frontal  gradients  which  are  about  an  order  of  magnitude  greater 
than  those  along  the  front.  Hoskins  (1974)  has  shown  that  similar  results  are  obtained  for 
the  steady  state  problem  as  well,  and  that  frontogenesis  increases  for  this  case  when  the 
thermal  wind  opposes  the  basic  flow. 

Although  the  Primitive  Equations  (PEs)  have  been  the  primary  tool  for  integrations 
used  in  numerical  weather  prediction,  they  have  the  disadvantage  of  being  so  general  as  to 
preclude  significant  analytical  manipulations.  In  particular,  the  lack  of  balance  conditions 
implies  the  desired  fields  (e.g.  wind  and  mass)  cannot  be  obtained  by  a  simple  invertibility 
relationship.  On  the  other  hand,  the  relatively  extreme  simplification  of  the  PEs  resulting 
in  the  Quasi-geostrophic  (QG)  approximation  cannot  mathematically  handle  the  dynamics 
which  occur,  among  other  areas,  in  and  near  a  midlatitude  front.  Fortunately,  the  study  of 
fronts  benefits  from  the  application  of  another  more  modern  tool  described  below  which 
has  been  successfully  applied  to,  for  example,  ageostrophic  secondary  frontal  flows  by 
Hoskins  and  Draghici  (1977)  and  the  balanced  squall  line  model  of  Schubert  et  al.  (1989). 
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The  Geostrophic  Momentum  Approximation  (first  alluded  to  by  Eliassen,  1943,  and 
hereafter  referred  to  as  GMA)  represents  a  compromise  between  the  complex  accuracy  of 
the  PEs  and  the  simplicity  of  QG  theory.  The  GMA  equations,  a  filtered  set  in  that  they 
cannot  describe  gravity  wave  propagation,  are  more  general  than  the  QG  equations  since 
fewer  terms  are  dr'pped  from  the  PEs.  In  fact,  they  bear  more  formal  resemblance  to 
the  PEs  than  they  do  to  those  resulting  from  QG  theory:  advection  is  by  the  total  wind 
and  vertical  advection  is  retained  in  the  momentum  equation.  Perhaps  most  importantly, 
though,  the  essential  properties  of  a  three  dimensional  vorticity  equation,  conservation  of 
potential  vorticity,  and  the  energy  equation  are  all  formally  consistent  with  those  derived 
from  the  PEs.  Additionally,  the  GMA  is  valid  for  flows  where  the  vertical  component  of 
relative  vorticity  is  on  the  order  of  the  coriolis  parameter — an  important  consideration  for 
the  study  of  fronts,  jets,  and  many  other  observable  atmospheric  phenomena.  It  should  be 
noted  here,  however,  that  the  GMA  is  not  valid  in  regions  where  the  flow  is  highly  curved 
(defined  as  any  region  where  the  curvature  vorticity  is  on  the  order  of  the  shear  vorticity) 
and  it  is  not  clear  what  if  any  effects  are  neglected  due  to  this  assumption. 

Ideally,  we  would  like  to  retain  the  benefits  of  a  general  theory  such  as  the  GMA,  while 
retaining  as  many  analytical  and  numerical  benefits  as  possible  arising  from  the  relative 
simplicity  of  QG  theory.  Hoskins  (1975)  and  Hoskins  and  Draghici  (1977)  made  progress 
toward  such  a  theory  by  advocating  a  coordinate  transformation  which,  when  applied  to 
the  GMA,  were  termed  the  “Semigeostrophic”  (SG)  equations.  With  these  new  horizontal 
coordinates,  the  advecting  velocities  become  geostrophic  as  in  the  QG  equations.  However, 
because  the  new  equations  are  mathematically  identical  to  the  GMA,  we  conclude  the 
ageostrophic  advections  are  retained  implicitly  in  the  coordinate  transformation.  Similarly, 
newly  defined  ageostrophic  components  allow  retention  of  both  a  continuity  equation  and 
vector  momentum  equation  which  are  formally  identical  to  QG  theory,  but  with  the  latter 
equation  retaining  vertical  motion  which  is  again  implicit  in  the  change  of  coordinates. 
The  thermodynamic  equation  becomes  more  general  since  its  vertical  motion  is  against 
the  (geostrophic)  potential  vorticity,  which  varies  in  three  dimensional  space,  instead  of 
being  against  the  static  stability  parameter,  which  varies  only  in  the  vertical. 
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Further  development  of  SG  theory  was  made  by  Schubert  et  al.  (1987,  1989)  with 
the  simultaneous  use  of  geostrophic  and  isentropic  coordinates,  which  allows  the  system  of 
equations  to  be  reduced  to  a  predictive  equation  for  the  potential  pseudo-density  (inverse 
Rossby-Ertel  potential  vorticity)  and  a  diagnostic  relationship  for  the  Bernoulli  function. 
With  this  choice  of  coordinates,  the  ageostrophic  circulation  becomes  entirely  implicit  as 
well.  Of  course,  for  adiabatic  motions  the  total  ^-derivative  disappears,  and  the  potential 
density  becomes  a  conservative  quantity.  The  potential  density  equation  has  only  one 
component  with  0  as  the  vertical  coordinate  since,  as  noted  by  Haynes  and  McIntyre 
(1987)  for  the  PEs  and  later  applied  to  the  GMA  equations  by  Schubert  et  al.  (1989), 
there  is  no  net  transport  of  potential  vorticity  across  any  isentrope;  thus,  potential  vorticity 
can  be  neither  created  nor  destroyed  in  a  layer  bounded  by  two  isentropes.  The  application 
to  the  GMA  equations  of  the  impermeability  of  isentropes  to  potential  vorticity  is  one  of 
several  indications  of  the  consistency  of  the  GMA. 

There  are  other  benefits  gained  by  utilization  of  the  SG  equations.  For  example,  be¬ 
cause  the  Jacobian  of  the  coordinate  transformation  is  also  the  dimensionless  vertical  com¬ 
ponent  of  absolute  vorticity,  resolution  of  a  numerical  scheme  is  enhanced  in  geostrophic 
space  in  regions  of  a  relative  maximum  of  this  quantity  (in  such  regions,  the  shift  obtained 
by  transforming  a  field  back  to  physical  space  is  the  greatest).  Thus,  for  a  given  grid  spac¬ 
ing,  features  which  may  be  unresolveable  in  physical  space  are  captured  in  geostrophic 
space  and  shifted  back  to  physical  space.  This  procedure  often  produces  realistic  replica¬ 
tions  (e.g.  vorticity  patterns  near  a  cold  front  produced  by  Hoskins  and  West,  1979). 

So  far  we  have  considered  from  a  historical  perspective  some  of  the  conditions  under 
which  upper  fronts  may  form,  and  discussed  the  relevance  of  certain  tools  and  approaches 
to  the  study  of  the  frontogenesis  process.  The  approach  of  Hoskins  mentioned  earlier 
has  aspects  which  preface  the  possibility  of  further  research.  Keyser  and  Pecnick  (1985) 
allude  to  these  possibilities  by  pointing  out  that  the  primary  success  of  two-dimensional 
semigeostrophic  theory  in  the  context  of  frontogenesis  has  thus  far  been  with  regard  to 
surface  fronts.  They  proceed  to  develop  a  PE  model  which  simultaneously  includes  the 
effects  of  confluence  (horizontal  deformation)  and  horizontal  shear,  along  with  three-  di¬ 
mensional  characteristics,  resulting  in  the  best  representations  of  upper  fronts  to  date  in  a 
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two-dimensional  model;  still,  there  is  the  lingering  question  of  whether  a  balanced  model 
can  successfully  model  internal  frontogenesis  in  two  dimensions  or,  as  stated  by  Buzzi  et 
al.  (1981)  with  regard  to  semigeostrophic  theory,  a  three  dimensional  representation  is 
required. 

Figure  1.1  shows  a  model-produced  upper  front  which  apparently  is  as  realistic  a 
result  as  has  been  achieved  for  two-dimensional  balanced  flows  to  date.  Notice  the  front 
doesn’t  extend  deep  into  the  troposphere,  and  the  slight  folding  of  the  dynamic  tropopause 
is  achieved  only  by  the  use  of  the  geostrophic  coordinate  in  the  model.  This  should  be 
compared  with  figure  1.2,  which  is  representative  of  the  growing  base  of  observational  data 
on  tropopause  folding  events.  As  is  frequently  the  case,  the  fronts  in  this  figure  formed 
within  a  developing  mid-latitude  cyclone,  with  a  fold  on  either  side  of  the  upper  trough  in 
the  height  field.  Here  we  see  the  formation  of  deep  fronts  with  well-defined  folds,  strong 
jet  streaks,  and  sharply  sloping  isentropes  in  the  vicinity  of  the  fronts. 

Thus  we  attempt  to  both  simplify  and  improve  on  previous  results  regarding  the 
frontogenesis  process.  Specifically,  we  address  the  following  questions:  (1)  Can  the  model 
equations  be  simplified  and  generalized  by  using  both  geostrophic  and  isentropic  coordi¬ 
nates?  (2)  Will  the  simultaneous  inclusion  of  (a)  a  continuous  potential  vorticity  field  in 
the  vertical  and  (b)  a  height-dependent  deformation  field,  which  has  not  been  attempted 
within  the  context  of  semigeostrophic  theory,  result  in  a  more  realistic  depiction  of  the 
upper- level  frontogenesis  process? 

In  order  to  answer  these  questions,  we  use  the  /-plane  system  of  equations  with 
the  GMA  and  proceed  in  Chapter  2  to  develop  the  three  dimensional  potential  density 
equation  and  associated  invertibility  relation  in  geostrophic  coordinates.  This  system  is 
closed  by  simply  including  the  geostrophic  and  hydrostatic  relations  which,  already,  is  a 
simplification  of  the  system  in  its  various  “physical”  coordinate  forms.  In  Chapter  3,  the 
reduction  of  this  model  to  two  dimensions  in  the  transformed  cross-frontal  plane  allows 
an  analytic  solution  of  the  prediction  equation  and  a  de-coupling  of  the  predictive  and 
diagnostic  equations  from  a  numerical  standpoint.  Frontogenesis  is  idealized  following  the 
general  approach  of  Hoskins  et  al.  (1971,  1972)  and  Buzzi  et  al.  (1981),  but  with  the 
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Figure  1.1:  A  two-dimensional  balanced  model  result  showing  formation  of  a  folded 
tropopause  or  constant  potential  vorticity  surface  (shown  as  the  heavy  solid  line)  and  asso¬ 
ciated  upper-front,  using  a  deformation  held  as  the  forcing  mechanism  and  the  geostrophic 
horizontal  coordinate.  Thin  solid  lines  are  isentropes,  and  isotachs  are  dashed.  Taken  from 
Hoskins  (1972). 


Figure  1.2:  Cross-section  through  an  Islandic  Low  as  observed  on  12  April  1983.  The 
tropopause  shown  is  a  constant  potential  vorticity  surface  and  is  represented  by  the  heavy 
solid  line  separating  stratospheric  air  (clear)  and  tropospheric  air  (shaded).  Thin  solid 
lines  and  dashed  lines  are  as  in  figure  1.1,  also.  Taken  from  Shapiro  (1985). 
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aforementioned  additions:  a  generalized  height  dependent  deformation  rate,  and  a  more 
realistic  initialization  of  the  potential  vorticity  field.  The  upper  and  lower  boundaries  are 
chosen  as  isentropic  as  well  as  isobaric  surfaces.  The  resulting  equations  are  shown  to 
take  the  same  form  as  those  for  the  specific  case  where  the  deformation  rate  is  constant. 
The  use  of  a  two  dimensional  model  has  the  historical  advantage  of  producing  realistic 
transverse  secondary  circulations  and,  in  addition,  with  our  choice  of  coordinates  we  show 
the  simplicity  with  which  these  circulations  can  be  recaptured. 

Results  show  the  combination  of  allowing  the  rate  of  deformation  to  vary  vertically 
and  prescribing  the  initial  field  of  potential  vorticity  as  a  continuous  one  do  indeed  produce 
realistic  upper  level  fronts,  perhaps  most  notably  a  true  folding  of  the  dynamic  tropopause 
which  in  some  cases  extends  deep  into  the  troposphere.  In  these  cases  the  tropopause  is 
defined  as  a  constant  potential  vorticity  surface  (see  Danielsen,  1980).  Although  at  this 
point  it  was  unclear  what  aspect(s)  of  our  model  produced  realistic  frontal  parameters, 
we  felt  the  key  ingredient  was  the  vertical  wind  shear  characteristic  of  the  deformation 
field.  Thus  an  additional  question  was  whether  a  generalization  of  the  wind  field,  wherein 
assumed  vertical  shears  are  consistent  with  baroclinic  wave  development,  would  produce 
similar  results.  Further  model  runs  in  Chapter  4  verify  this  hypothesis  with  respect  to 
formation  of  upper  fronts.  In  this  case  we  also  redefined  the  initial  state  in  a  manner 
known  to  isolate  formation  of  internal  fronts;  thus,  results  of  Chapter  4  verify  that  realistic 
modeling  of  upper  frontal  parameters  can  be  accomplished  when  disassociated  from  surface 
fronts. 

Chapter  2  reviews  the  three-dimensional  theory  from  which  our  model  is  derived. 
Chapter  3  defines  the  two-dimensional  model,  describes  results  of  selected  model  runs,  and 
compares  these  results  with  previous  work.  In  Chapter  4  further  testing  and  comparison 
is  done  by  redefining  the  initial  field.  A  summary  follows  in  Chapter  5. 


Chapter  2 


THEORY  REVIEW 

2.1  Review  of  the  Potential  Density  Equation 

We  begin  with  the  /-plane  system  of  equations  with  the  geostrophic  momentum 
approximation  (Eliassen,  1948;  Hoskins,  1975).  Proceeding  with  an  analysis  similar  to 
that  of  Schubert  et  al.  (1989),  we  will  arrive  at  a  simple  predictive  equation  for  the 
potential  pseudo-density  (inverse  Rossby-Ertel  potential  vorticity).  Later,  in  our  two 
dimensional  model,  this  predictive  equation  will  take  on  an  even  simpler  form.  Assuming 
the  flow  is  frictionless,  and  using  potential  temperature  as  the  vertical  coordinate,  our 
system  becomes 
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where  ( ug,va )  =  f~l(-dMJOy,dM./dx)  are  the  components  of  geostrophic  velocity,  (u,u) 


the  horizontal  components  of  the  total  velocity,  II  =  cp(p/po)K  the  Exner  function,  M  = 
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theta  velocity,  and 

n  a  A  a  .  A 

(2.5) 


D  d  d  d  id 
Dt  dt  T  Ox  r  dy  09 


the  total  derivative.  The  vertical  component  of  the  vorticity  equation  can  be  derived  from 
(2.1)  and  (2.2)  and  takes  the  form 
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where  the  components  of  the  absolute  vorticity  vector  are  given  by 
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and  X  and  Y  are  defined  below  in  (2.10).  We  can  eliminate  the  divergence  between  (2.4) 
and  (2.6)  to  obtain  one  form  of  the  potential  density  equation: 
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where 

a*  =  Lj  (2.9) 

is  the  potential  pseudo-density,  or  potential  density  (PD)  for  short.  With  9  as  the  vertical 
coordinate  we  can  simplify  our  system  of  equations  (hence,  the  vorticity  and  pseudo- 
density  equations  also)  significantly  for  motions  which  are  adiabatic.  We  will  make  the 
adiabatic  assumption  later  in  chapter  3  when  we  develop  the  (dry)  frontal  model. 


2.2  Review  of  Semigeostrophic  Theory 
2.2.1  Coordinate  Transformation 


Since  (2.8)  involves  advection  by  the  total  wind,  this  equation  along  with  the  invert- 
ibility  relation  we  formulate  in  section  2.2.2  will  not  form  a  closed  system.  Expressed  in 
physical  space,  then,  this  system  would  not  be  appropriate  for  modeling  purposes  unless 
we  were  to  solve  an  additional  equation.  With  this  as  the  initial  motivation,  we  now  intro¬ 
duce  the  geostrophic  coordinates,  where  in  this  case  we  simultaneously  use  the  isentropic 
vertical  coordinate  approach  following  the  arguments  of  Schubert  et  at.  (1989).  The  hori¬ 
zontal  transformation  has  the  affect  of  switching  the  function  of  physical  and  geostrophic 
coordinates,  respectively,  by  making  geostrophic  variables  independent  variables  and  the 
physical  variables  dependent  variables.  We  define  geostrophic  coordinates  in  the  usual 
way: 

(X,  Y,  0,T)  =  (x  +  vgff,  y  -  us//,0,  t),  (2.10) 

where  the  term  “geostrophic”  is  a  result  of  the  fact  that,  for  the  frictionless  case  described 
in  our  governing  equations,  DX/D t  =  ug  and  DY/D t  =  vg,  as  can  be  easily  seen  from 


s 


11 

(2.1),  (2.2),  and  (2.10).  Thus,  a  particle  following  the  motion  in  geostrophic  space  moves 
with  the  geostrophic  wind. 

The  relationships  between  derivatives  in  (x,y,0)  and  (X,  Y,  0)  space  are  written  as 


d 

dX  d  dY  d  d 

(2.11) 

dt 

dt  dX  +  dt  dY  +  3T’ 

d 

dX  d  dY  d 

(2.12) 

dx 

dx  dX  +  dx  dY  ' 

d 

dX  d  dY  d 

(2.13) 

dy 

dy  dX  +  dy  dY  ’ 

d 

dX  d  dY  d  d 

(2.14) 

dd  ~ 

dO  dX  +  dO  dY  +  30* 

The  transformations  (2.11)-(2.14)  imply  the  total  derivative  operator  (2.5)  can  be  written 

(2.15) 
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Thus,  (2.8)  is  still  the  fundamental  predictive  equation  of  the  model  but  with  the  total 
derivative  now  given  by  (2.15).  We  can  invert  (2.12)  and  (2.13),  which  yields 


C  9 

dY  d  dY  d 

(2.16) 

fdX 

dy  dx  dx  dy' 

C  d 

dXd  dXd_ 

(2.17) 

fdY 

dy  dx  dx  dy' 

and 


We  can  use  (2.16)  and  (2.17)  in  (2.14),  along  with  the  Jacobian  relations  (2.7),  to  obtain 

(2.18) 
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We  see  from  (2.18)  that  the  term  “vortex  coordinates”  can  reasonably  be  applied  to  the 
geostrophic  coordinates,  since  the  partial  theta-derivative  in  transformed  space  is  along 
the  three  dimensional  absolute  vorticity  vector  in  ( x,y,0 )  space. 

An  important  aspect  of  the  geostrophic  momentum  approximation  is  that  the  essential 
conservation  relations  of  the  primitive  equations  are  maintained  in  their  original  form. 
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For  example,  conservation  of  potential  temperature  is  formally  identical  to  the  primitive 
equation  form,  but  with  the  total  derivative  given  by  (2.15);  a  similar  statement  holds  for 
the  energy  equation  as  well. 

It  is  easy  to  imagine  from  the  form  of  the  potential  density  equation  that  its  inverse, 
the  potential  vorticity  equation,  will  be  formally  identical  to  the  primitive  equation  form 
of  the  potential  vorticity  equation.  In  fact,  the  potential  vorticity  equation  is  derived  in 
the  same  manner  as  the  potential  density  equation  (2.8).  Thus,  combining  (2.4)  and  (2.6) 


can  also  yield 


Df  a  \  dx 


(2.19) 


where  q  =  £/<r  is  the  Rossby-Ertel  Potential  Vorticity.  Notice  that  in  (x,y,9),  only  the 
vertical  component  of  vorticity  contributes  to  the  potential  vorticity1  equation  or,  in  the 
case  of  (2.8),  the  potential  density  equation  .  The  reason  for  this  is  related  to  the  fact 
that  (2.6)  can  be  written  in  the  form 


d(aq)  d(uoq-6 d  (yog  -  07?) 


so  that  intuitively  we  see  that  what  Haynes  and  McIntyre  (1987)  proved  for  the  isentropic 
form  of  the  primitive  equations  holds  as  well  when  the  geostrophic  momentum  approx¬ 
imation  is  made;  namely,  isentropes  are  semi-impermeable  to  potential  vorticity  (hence, 
to  potential  density  also)  in  the  sense  that  there  can  be  no  net  transport  of  q  across  any 
isentrope.  Applying  (2.18)  to  (2.19)  yields 


D  q  C  00 
Df  q  d 0 


(2.20) 


which  is  again  considerably  simpler  than  its  counterpart  in  ( x,y,6 ).  Similarly,  we  can 
write  (2.8)  in  the  form 

5£l  +  <T*J^=0.  (2.21) 

D  t  de 


'This  vorticity  quantity  is  more  properly  referred  to  as  the  isentropic  vorticity  (g,  to  be  distinguished 
from  the  vertical  vorticity  component  in  pressure  or  height  coordinates.  However,  we  will  keep  the  con¬ 
vention  of  using  the  symbol  (  when  referring  to  the  vertical  vorticity  component  in  all  systems. 
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Already,  important  features  of  geostrophic  coordinates  are  apparent.  In  addition  to  the 
simplicity  of,  for  example,  (2.20),  the  Jacobian  of  the  transformation  from  physical  co¬ 
ordinates  is  simply  the  dimensionless  vertical  component  of  absolute  vorticity  £//  (com¬ 
pare  (2.12)  and  (2.13)  with  (2.7)).  Second,  horizontal  advecting  velocities  have  become 
geostrophic,  which  eliminates  the  need  to  solve  an  additional  equation  to  close  the  sys¬ 
tem  in  the  adiabatic  case,  yet  the  governing  equations  implicity  retain  the  ageostrophic 
advections.  In  our  two-dimensional  frontogenesis  model,  a *  will  be  the  predicted  quantity. 

2.2.2  Invertibility  Principle 


Because  prediction  of  a*  can  now  be  performed  in  geostrophic  coordinate  space,  in 
order  to  form  a  closed  system  we  need  to  follow  the  procedure  of  Schubert  et  al.  (1989) 
and  formulate  the  invertibility  relation  alluded  to  above  in  transformed  space  also.  It  can 
be  shown  that  the  geostrophic  and  hydrostatic  relations  in  (X,  Y,  0)  take  the  form 


dU*  0M*  dU1 
dX  '  dY  '  <90 


)• 


(2.22) 


which  is  identical  to  the  form  taken  in  ( x,y,$ )  space  with  M  replaced  by  M*,  where 


M* 


(2.23) 


is  the  Bernoulli  function.  The  potential  density  <7*  is  a  combination  of  the  mass  field  a  and 
the  wind  field  £.  However,  a  is  related  to  II,  and  using  (2.7)  C  can  be  expressed  in  terms 
of  ug  and  vg.  By  (2.22),  then,  the  wind  and  mass  fields  can  be  obtained  from  <7*  if  we 
can  somehow  invert  it  to  obtain  M*.  This  second  order  partial  differential  equation,  along 
with  appropriate  boundary  conditions,  is  usually  referred  to  as  the  invertibility  principle. 
A  relation  between  M*  and  <r*  can  be  derived  as  follows.  From  the  definition  of  <r*  we 
have 

!f£  +  r<r*  =  °,  (2.24) 

where  T  =  dU/dp  =  ull/p.  This  last  equation  can  be  written 
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Expressing  x  and  y  in  terms  of  ug  and  vg  by  (2.10),  and  then  using  (2.22),  we  can  write 
(2.25)  as 

f2  d2Mm\fd2Mm\2  /2  d2M*\/S2M*\2  d2M'\  /  d2M*\  f  82W\ 

y  dx2  J  \dYd9 J  +  y  dY2  )  { dxdo)  +2[dxdYj  \dYdo)  [ aeax j 

\  2j  2  d2W  d2M*\  (d2W\(d2U'\  /a2M*\2la2M* 

+  y  dx 2  dY2  )  +  {  ax2  J  [d  Y2  J  (dXdY )  J  dQ2 

+/4rcr*  =  0.  (2.26a) 

If  the  upper  boundary  is  an  isentropic  and  isobaric  surface  with  potential  temperature  ©t 
and  Exner  function  IIt,  the  upper  boundary  condition  for  (2.26a)  is  simply 

5M* 

~~  =  nT  at  0  =  0T.  (2.26b) 


If  we  neglect  the  effects  of  topography  and  assume  that  the  lower  boundary  is  both  the 
constant  isobaric  height  surface  0  =  0  and  the  isentropic  surface  0  =  0b,  then  M  =  011 
at  0  =  0g.  Written  in  terms  of  M*,  this  lower  boundary  condition  becomes 


0  =  0B- 


(2.26c) 


The  lateral  boundary  conditions  depend  on  the  particular  application  and  will  be  discussed 
later  in  the  context  of  the  two  dimensional  model.  It  should  be  noted  here  that  defining  the 
lower  boundary  as  both  isobaric  and  isentropic  does  not  necessarily  exclude  formation  of  a 
surface  front;  however,  in  Chapter  4  we  prescribe  these  boundaries  to  be  constant  surfaces 
of  potential  density;  in  this  case,  as  noted  by  Keyser  et  al.  (1986),  the  possible  coupling 
effects  between  surface  and  upper-level  fronts  (through  vertical  motions)  are  excluded. 

Equations  (2.8),  (2.22)  and  (2.26)  form  a  closed  system.  The  computational  scheme 
is  as  follows:  knowing  <?*,  solve  (2.26)  for  M*;  use  (2.22)  to  compute  ug  and  vg ;  use  these 
geostrophic  winds  in  (2.8)  to  predict  a  new  a*. 


Chapter  3 


FRONTOGENESIS  BY  HEIGHT  DEPENDENT  DEFORMATION 

3.1  Two  Dimensional  Model 

3.1.1  Theory 

So  far,  using  both  the  geostrophic  momentum  approximation  and  a  useful  coordi¬ 
nate  transform,  we  have  reviewed  the  necessary  three  dimensional  theory  to  study  a  wide 
range  of  fluid  phenomena  which  are  too  general  to  be  studied  with  quasigeostrophic  the¬ 
ory;  however,  we  have  at  the  same  time  retained  its  major  benefits,  such  as  the  balance 
relationships  in  the  resulting  equations  which  are  necessary  to  recover  the  fields  of  inter¬ 
est.  We  now  apply  a  simplified  version  of  this  theory  to  study  upper  level  (i.e.  internal) 
frontogenesis.  In  particular,  we  assume  the  motions  to  be  adiabatic,  and  we  simplify  to 
two  dimensions  both  the  potential  density  equation  and  the  invertibility  relation  arrived 
at  in  Chapter  2  in  order  to  study  the  frontogenesis  process  in  the  transverse  (X,0)  plane. 
The  assumption  that  no  variations  in  the  fields  of  interest  occur  along  the  front  follows 
a  frequently  held  axiom  that  the  relationship  between  such  variations  and  the  frontogen¬ 
esis  process  itself  is  not  clear,  but  that  such  an  assumption  is  acceptable  for  along-front 
geostrophic  flows  (i.e.  transverse  geostrophy)  where  the  front  is  assumed  to  be  straight 
and  variations  of  the  along-front  fields  are  small. 

Let  us  now  generalize  the  two-dimensional  frontogenesis  problem  of  Hoskins  (1971, 
1972)  and  Hoskins  and  Bretherton  (1972).  In  this  case,  fronts  oriented  in  the  y-direction 
are  assumed  to  be  forced  by  a  pure,  height  dependent  deformation  field  so  that 


ug(x,0)  =  -a(0)x, 

(3.1) 

vg(x,y,0,t)  =  a{0)y  +  v'g{x,0,t). 

(3.2) 
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where  a  =  a(0)  is  the  vertically  varying  deformation  rate,  the  first  terms  on  the  right 
hand  side  represent  the  fixed  (or  “slowly”  varying)  deformation  field,  and  the  v'  term 
represents  the  rotational  flow  generated  during  frontogenesis  so  that  i^(x,0, 0)  =  0. 

Since  we  assume  that  motions  along  the  front  are  geostrophic,  the  entire  secondary 
circulation  is  in  the  transverse  plane  and  is  given  by  ua.  We  can  easily  recover  this 
circulation;  from  the  definition  u  —  Dx/D t,  equations  (2.10),  the  geostrophic  relation 
(2.22),  the  assumption  (3.1),  and  the  fact  that  ua  =  u  -  ug,  we  have 


ua 


P 


d  dM *  v02M’ 
ST  dX  ~  aX  dX 2 


(3.3) 


3.1.2  Governing  Equations 


Now  we  formulate  the  two  dimensional  invertibility  principle  in  a  manner  consistent 
with  the  requirements  of  the  numerical  solver  which  is  discussed  later.  First,  however,  we 
reduce  our  fundamental  predictive  equation  (2.8)  to  two  dimensions  in  geostrophic  and 
isentropic  space.  Setting 

d_  _  dXd_  dYJP 
dy  dy  dX  dy  dY 

with  respect  to  the  field  variables,  using  the  definitions  (2.10)  and  the  assumptions  (3.1) 
and  (3.2),  we  obtain 

JL  =  _“JL  (3  4) 

dY  fdx'  [  ; 

Using  the  adiabatic  assumption  0  =  0  and  substituting  (3.4)  into  (2.15)  transforms  our 
total  derivative  operator  into 

n  f)  A 

(3.5) 


d.  =  A_aXA 

d*  dt  ex' 


so  that  equation  (2.8)  becomes 


OT  X  dX  °’ 


which  has  as  its  solution 


a'  (X,0,T)  =  af  (XeoT,e)  . 


(3.6) 


(3.7) 
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Substituting  (3.4)  into  (2.26a),  we  obtain  the  corresponding  two-dimensional  invertibility 
relation 


f  d2M*  \ 
(^0dX  ) 


d2M*  d2M* 
dX2  <90 2 


+  /2 


=  0. 


1  lUCOC  U.C11L  l/i 


Since  a  is  coupled  to  T  in  the  prediction  of  o’,  only  the  product 
important;  in  particular,  we  can  make  a  as  small  as  is  convenient  for  our  purposes.  If  we 
choose  a  so  that  a2//2  <C  1,  then  we  can  approximate  the  above  as 


d2M* 

<902 


(/2-^)  +  (0x)2  +  /^‘=°-  <3-8a> 


which,  interestingly,  is  identical  to  the  form  obtained  for  both  the  specific  case  where  a  is 
considered  constant,  and  for  the  case  study  of  the  atmospheric  response  to  a  moving  squall 
line  (Schubert  et  al.  1989)  where  both  the  zonal  geostrophic  and  meridional  ageostrophic 
components  of  the  wind  are  assumed  to  be  zero.  The  upper  boundary  condition  is  un¬ 
changed: 

=  nT  at  0  =  0T,  (3.8b) 


and  from  (2.26c)  and  (3.4)  the  lower  boundary  condition  becomes 

...  1  fd  n 

0  80"'“  +27i(«X  J  "°  "  0-6b’ 


(3.8c) 


where  we  again  use  the  constraint  a2//2  <C  1.  For  the  lateral  boundary  condition  we 
make  use  of  the  fact  that  the  geostrophic  deformation  field  (3.1)  and  (3.2)  implies  vg  — *• 
OasX  -*■  dboo.  In  the  model,  this  last  expression  will  be  approximated  by 


M*  (0,T)  =  M*  (0,0)  at  X  =  -Xl,3Xl,  (3.8d) 

where  we  use  the  definition  (2.22)  and  define  Xl  as  a  (variable)  scale  width  of  the  domain 
of  interest  (see  the  numerical  procedures  section). 

From  a  practical  standpoint,  we  can  see  that  already  our  two  dimensional  model  is 
much  simpler  to  solve  and  requires  less  computer  time  due  to  the  analytic  form  of  the 
solution  (3.7).  The  advecting  velocity  in  (3.6)  is  simply  the  product  of  the  independent 
variable  X  and  the  deformation  rate  a,  which  is  prescribed  as  varying  only  in  the  vertical; 
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therefore,  we  need  not  diagnose  M*  in  order  to  determine  this  advection  quantity  as  we  did 
for  the  three  dimensional  case.  Thus,  in  two  dimensional  (X,  0)  space  our  computational 
scheme  is  as  follows:  knowing  a*,  use  (3.8)  to  diagnose  M*  only  for  those  time  steps 
for  which  we  wish  to  determine  the  wind  and  mass  fields.  To  recover  the  transverse 
ageostrophic  circulation  requires  diagnosing  M*  for  three  successive  time  steps  in  order  to 
make  use  of  a  discretized  version  of  (3.3). 

We  should  point  out  that  an  important  limitation  of  our  two-dimensional  model 
results  by  neglecting  along- front  temperature  gradients.  Keyser  and  Pecnick  (1985)  have 
shown  these  gradients  to  be  significant  in  the  formation  of  some  upper  fronts.  From  this 
point  of  view,  then,  our  results  should  be  taken  within  the  constraints  of  the  assumptions 
made  when  formulating  the  simplified  model  presented  here. 

3.2  Procedure  for  the  Invertibility  Solver 


Our  two  dimensional  model  has  been  reduced  to  (3.6)  and  (3.8),  the  former  having  an 
analytical  solution  given  by  (3.7).  The  numerical  solution  to  the  nondimensional  form  of 
(3.8)  developed  below  will  be  given  by  the  Full  Multigrid  (FMG)  solver  of  Fulton  (1989), 
an  efficient  program  utilizing  nonlinear  Gauss-Seidel  relaxation  for  smoothing.  The  FMG 
program  requires  that  Af(X,0)  =  M*(X,  0)  -  M  (0)  be  a  specified  value  on  each  lateral 
boundary;  that  is,  the  model  requires  Dirichlet  lateral  boundary  conditions1.  For  a  more 
in-depth  discussion  of  the  invertibility  solver  the  reader  is  referred  to  Fulton  (1989). 

To  rewrite  our  diagnostic  equation  (3.8)  in  the  nondimensional  form  needed,  we  in¬ 
troduce  the  nondimensional  coordinates  (it  should  be  noted  that  these  are  quite  different 
from  the  physical  coordinates  x  and  z) 


x  ~ 


0-0B 
0T  ~  ©B  ’ 


(3.9) 


I 


'Here,  M*  is  a  horizontally  averaged  value  of  M*  and  At  is  a  deviation  value.  To  determine  an 
appropriate  value  for  M  ,  note  that  by  integrating  (3.13)  we  can  write  a  basic  state  form  of  the  initial 
pressure  field.  Integration  of  the  hydrostatic  relation  (2.22)  and  the  definition  of  II,  then,  will  define  a 
value  of  M  for  any  point. 
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>» 


with  c  =  [<zR(0t  —  ©b)]1^2,  <*  =  (pB  —  Pt)/pb,  and  R  the  gas  constant  for  dry  air. 
Derivatives  between  these  nondimensional  quantities  and  those  in  geostrophic  space  are 
then  related  by 


3 

dX 


LJL 

c  dx ' 


d 

dQ 


d 


(3.10) 


(0T  -  0b)  dz 

Now  we  write  (3.8)  in  its  nondimensional  form:  nondimensionalizing  M  by  c2,  a *  by  cr0 
(defined  in  the  next  section),  T  by  To  =  R/pb>  n  by  cp,  and  p  by  pe,  and  retaining  the 
original  notation  for  these  quantities  yields 


1  - 


d2M\  /=_  d2Ms 
dx 2  J  \Ta  dz 2 , 

l-(£) 

(3.11a) 

dM 

dz  at 

2  =  1, 

(3.11b) 

dM  1  SdM\2 

dz  2\dx) 

=  0  at  2  =  0, 

(3.11c) 

M  =  0  at  X 

=  -Xl  ,3Xl  , 

(3. lid) 

where  b  ~  ©b/(©t  -  ©b).  Second-order  centered  finite  differencing  is  used  for  both 
(3.11a)  (the  domain  interior),  and  for  (3.11b)  and  (3.11c)  (the  upper  and  lower  boundary 
conditions)  where  ghost  points  satisfy  the  requirements  of  the  finite  difference  scheme. 
Applying  the  homogeneous  characteristics  of  the  atmosphere  outside  the  domain  allows 
use  of  Dirichlet  lateral  conditions  (discussed  below)  for  the  four  ghost  points  on  those 
boundaries.  Once  the  Ad-field  has  been  diagnosed,  M*  may  be  easily  found  and,  using 
(2.22),  the  wind  and  pressure  fields  calculated.  Representation  of  the  fields  of  interest 
in  physical  space  can  be  made  by  using  (2.10)  once  the  geostrophic  wind  field  has  been 
calculated. 

If  there  are  nz  vertical  grid  intervals,  the  vertical  resolution  is  A0  =  (0t  -  ©b  )/nz 
K.  To  determine  the  horizontal  grid  spacing,  we  first  must  examine  the  horizontal  extent 
necessary  to  avoid  contamination  of  the  results  by  the  lateral  boundary  conditions,  defined 
as  one  Rossby  length  A,  where 


A  = 


©T  -  ©B 

/ 


(R/pb)(pb  -  Pt) 


©T  -  ©B 


(3.12) 


s 
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From  this  information  we  define  Xl  =  2A,  so  that  the  domain  is  eight  Rossby  lengths 
wide.  Although  this  domain  size  is  somewhat  cumbersome  from  a  numerical  standpoint, 
the  realism  of  our  results  here  demand  that  warm  air  border  the  cold  trough  on  both  sides. 
The  asymmetry  of  the  domain  definition  allows  for  a  shifted  deformation  axis  as  described 
below.  With  our  values  for  pressure  and  potential  temperature  given  in  the  next  section 
and  /  taken  as  10-4  s_1,  A=1595  km;  the  relatively  large  value  for  A  is  a  result  of  the 
depth  of  the  domain.  The  horizontal  resolution  is  given  by  AX  =  A[A0/(0x  —  0b)]-  A 
typical  value  for  nz  is  64,  which  corresponds  to  A0=1.56  K  and  AX=25  km. 

3.3  Description  of  the  Initial  Domain 


Equations  (3.6)  and  (3.8)  together  with  the  geostrophic/hydrostatic  relations  (2.22) 
form  a  closed  system;  However,  diagnosis  of  the  wind  and  mass  fields  by  (3.8)  requires 
that  the  initial  potential  density  field  be  prescribed  in  (X,0)  space.  To  gain  some  insight 
as  to  the  observed  vertical  variation  of  potential  density,  we  derived  the  a *  field  using  data 
from  the  U.S.  Standard  Atmosphere  tables  (1966)  for  the  months  of  January  and  July. 
Based  on  this  data,  we  define  the  initial  pressure  field  as 


Pl(X,0)  =  pt  + 


(<J^)(PB~PT)  +  (*0HlftA 
PB  -  PT  +  ( SQScr )  In  B 


(pb  -  Pt)  , 


for  which  the  initial  potential  density  field  2  is 


*i(X,0) 


<7o  -  (tfcr)tanh 

1  +  f  J&sL)  In  B 

\PB~PT  / 


(3.13) 


(3.14) 


2The  relationship  between  the  initial  fields  of  pseudo-density  and  potential  pseudo-density  is  found  as 
follows:  with  the  geostrophic  components  as  defined  above  in  (3.1)  and  (3.2),  we  know  initially  that 


dug 

dx 


=  -a, 


dvg 

and  —  =  o, 

dy 


so  that  using  (2.7)  we  can  write  Cl  =  /  —  °J//  =  /•  Thus,  from  the  definition  <r*  =  (// Qo  we  can  say 
with  a  high  degree  of  accuracy  that  initially  the  potential  density  field  and  the  density  field  are  equivalent 
in  the  domain,  and  the  two  quantities  are  essentially  interchangeable.  In  addition,  the  condition  (3.8d) 
implies  that  in  the  far  field  and  for  all  time  M*  F(X),  so  we  have  exactly  (  =  /  and  o’  =  a  outside  the 
domain. 
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where 

c°Sh  (g^1) ' 
cosh(^)J’ 

0B  =  280  K,  0T  =  380  K,  pB  =  101575  Pa,  pT  =  11575  Pa,  a0  =  (pB  -  Pt)/(0t  -  0B)  = 
900  Pa/K,  S<7  =  800  Pa/K,  and 

'  0  -  (A0)  [l  -  cos  j  -XL  <  X  <  3XL, 

01  =  _  (3.15) 

k  0  X  <  -XL  or  X  >  3XL, 

is  the  variation  of  potential  temperature  along  the  tropopause  in  the  (X,0)  plane.  The 

quantity  0  =  340  I<  is  a  global  average  value  of  the  tropopause  potential  temperature3.  As 

can  be  seen  in  figures  3.1  and  3.2,  with  the  above  definitions  the  initial  potential  density 

field  o\  closely  follows  the  observed  annual  average  (an  average  of  the  January  and  July 

values)  of  a.  Although  a  more  exact  fitting  to  observations  can  be  obtained,  the  values 

chosen  here  display  the  characteristic  sharp  decline  of  a  observed  in  regions  favorable  for 

frontogenesis — features  which  are  smoothed  out  in  the  climatological  averages  depicted  by 

the  observations. 

Thus,  unlike  the  two  layer  model  of  Hoskins  (1971),  we  prescribe  a  continuous  poten¬ 
tial  vorticity  (PV)  field  with  a  sharp  vertical  gradient  in  the  vicinity  of  the  tropopause. 
Notice  that  at  the  upper  and  lower  boundaries  (0t  and  0b),  equation  (3.13)  reduces  to 
the  constant  pressure  surfaces  px  and  pB  which  satisfy  our  original  boundary  assumptions 
in  equations  (2.26b)  and  (2.26c).  These  lower  boundary  values  were  chosen  as  an  annual 
average  from  the  climatological  data,  and  upper  boundary  values  were  chosen  in  such  a 
way  as  to  satisfy  our  desired  value  of  <7o-  The  parameter  A0  in  equation  (3.15)  controls 
the  initial  horizontal  temperature  gradient  in  the  vicinity  of  the  tropopause.  In  general 
terms,  its  numerical  value  is  small  in  keeping  with  the  approach  of  Hoskins  that  the  as¬ 
sumed  geostrophic  deformation  field  is  acting  on  a  thermal  gradient  which  is  initially  weak. 

3Referring  to  figures  3.1  and  3.2,  our  assumed  value  compares  acceptably  well  with  observed  tropopause 
temperatures.  Thu8,_for  the  purposes  of  our  model  we  consider  the  tropopause  to  have  a  globally  averaged 
value  equivalent  to  0,  and  within  the  domain  of  interest  the  tropopause  lowers  in  conjunction  with  the 
initial  horizontal  thermal  gradient  in  the  interior. 


cosh 

[c°5h  (fi^Si)J’ 


SUE 


model 


observed 


Figure  3.1:  Potential  density  vs.  potential  temperature.  Observed  values  are  for  the 
month  of  January  as  derived  from  the  U.S.  Standard  Atmosphere  tables.  The  dashed 
curve  is  the  representation  of  equation  (3.14)  for  a  typical  tropopause  temperature.  Un.  "3 
on  the  horizontal  axis  are  Pa/K. 
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Thus,  the  tropopause  has  a  trough  representing  colder  air  in  the  center  of  the  domain  in 
a  manner  similar  to  that  of  Hoskins  et  al.  (1985). 

Figures  3.3-3. 6  show  the  initial  state  of  the  entire  domain.  In  figure  3.3  the  a*  field 
is  greatest  at  the  surface  and  in  the  center  of  the  domain  below  the  upper  trough.  The 
associated  potential  vorticity  field  q  is  depicted  in  figure  3.4.  From  the  definitions  in 
Chapter  2,  we  have  the  relationship  q  =  //c*;  since  these  two  quantities  are  inversely 
proportional  within  the  constant  /,  emphasis  will  be  placed  on  the  more  familiar  q-field. 

Figure  3.5  confirms  our  earlier  result  with  regard  to  the  initial  vorticity  field;  the 
initial  meridional  geostrophic  wind  field  displayed  in  the  figure  indicates  a  very  weak 
horizontal  gradient  so  that,  as  stated  earlier,  the  earth’s  vorticity  is  nearly  equal  to  the 
absolute  vorticity.  Finally,  figure  3.6  gives  the  temperature  profile. 

3.4  Results  of  the  Model 

3.4.1  Initial  Results 

Preliminary  experiments  conducted  using  a  deformation  rate  which  linearly  decreased 
with  height  were  unsuccessful  at  producing  upper  fronts.  Desired  major  characteristics, 
most  importantly  a  folding  event,  were  not  reproduced.  Also,  neither  a  linear  nor  Gaussian 
definition  of  a  were  successful  when  the  zero  zonal  wind  line  (center  of  the  deformation 
field)  was  placed  in  the  center  of  the  trough.  As  noted  below,  however,  there  was  a  marked 
improvement  when  the  zero  line  was  moved  off  to  either  the  left  or  right  midway  between 
the  cold  and  warm  air.  Observed  features  of  baroclinic  waves  are  consistent  with  this  shift 
and  serve  as  sufficient  motivation;  the  existence  of  confluence  in  the  entrance  region  of  jet 
streaks  is  well  documented  both  upstream  and  downstream  of  the  associated  upper  trough, 
corresponding  to  a  model  shift  left  and  right,  respectively.  It  is  in  such  regions  that  the 
bulk  of  observational  data  has  been  taken  verifying  upper  frontogenesis  and  associated 
folding  events. 

3.4.2  Gaussian  Deformation 


Here  the  deformation  field  is  shifted  off  the  center  of  the  domain  to  a  point  halfway 
between  the  coldest  and  warmest  air.  There  are  two  such  points:  one  corresponds  to 
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Figure  3.6:  Initial  potential  temperature  field  in  (x,p)  space.  Isentropes  are  every  5  I<, 
ranging  from  the  lower  boundary  value  of  280  K  to  380  K  on  the  upper  boundary.  Distance 
is  given  in  km  on  the  horizontal  axis. 
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X  =  0,  and  the  other  to  X  =  2Xl-  The  shift  provided  by  far  the  most  realistic  results  of 
a  folded  dynamic  tropopause. 

In  this  experiment  we  set  a  =  10-sexp{-(l/2)[(©  -  0)/60]2},  where  0  is  the  level 
of  maximum  deformation  and  £0  is  the  e-folding  half-depth  of  the  a  field.  Typical  values 
for  these  quantities  were  330  K  and  20  K,  respectively.  The  coefficient  10-5  was  chosen 
so  that  our  assumption  (a//)2  <  1,  which  was  used  to  arrive  at  the  two-dimensional 
diagnostic  equation,  would  be  satisfied.  The  only  consequence  of  such  a  procedure  is  more 
time  must  elapse  to  achieve  the  same  results  than  would  have  been  the  case  for  a  larger 
value  of  a.  Although  references  to  time  can  therefore  only  be  made  in  a  relative  sense,  we 
nevertheless  prefer  it  to  references  to  the  product  qT,  which  is  somewhat  more  ambiguous 
since  a  varies  with  height. 

Figure  3.7  shows  that  the  folding  process  has  already  begun  in  the  PV  field  at  the 
24  hour  point.  A  jet  along  the  front  has  rapidly  developed,  as  shown  in  fig.  3.8  with 
maximum  winds  of  about  ^4  m/s  in  the  core  near  the  tropopause.  At  the  36  hour  point 
the  folding  process  has  intensified  with  a  deeper  region  of  anomalously  high  PV  values 
contained  in  the  developing  upper  front.  The  associated  jet  streak  has  increased  slightly 
in  magnitude  (figs.  3.9,  3.10).  After  48  hours  a  well  developed  folding  of  the  PV  field 
has  occurred;  figs.  3.11  and  3.12  show  the  PV  field  in  geostrophic  and  physical  space, 
respectively,  at  this  time.  There  are  several  noteworthy  features  in  comparing  the  two 
figures.  First,  the  region  above  the  fold  and  in  the  vicinity  of  the  tropopause  is  shown 
to  have  a  tighter  gradient  in  physical  space,  which  more  closely  matches  observations  of 
a  near  discontinuity  of  the  PV  field  near  the  tropopause.  Below,  the  fold  itself  displays 
the  same  characteristic,  and  is  also  more  narrow.  In  addition,  the  shift  back  to  physical 
space  has  increased  the  horizontal  extent  of  the  fold  by  moving  points  near  the  jet  farther 
west  than  corresponding  points  lower  in  the  troposphere.  All  of  these  features  are  a  result 
of  the  horizontal  and  vertical  wind  shear  which,  of  course,  is  greatest  near  the  jet.  Figs. 
3.13  and  3.14  show  the  wind  field  at  48  hours,  also  in  geostrophic  and  physical  space.  A 
more  realistic  (stronger)  wind  shear  pattern  results  from  the  shift.  In  general,  the  entire 
region  of  frontogenesis  has  been  tightened  due  to  its  cyclonic  nature.  The  characteristic  of 
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Figure  3.8:  Meridional  geostrophic  wind  field  at  24  hours  in  (X,0)  space.  Strongest, 
northerly  winds  are  12  m/s,  with  a  southerly  jet  core  of  44  m/s  (isotachs  every  4  m/s). 
Distance  on  the  horizontal  axis  is  in  km. 
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shrinking  low  pressure  regions  and  broadening  high  pressure  (anticyclonic)  regions  has,  as 
stated  previously,  been  well  documented  by  numerous  authors  of  semigeostrophic  theory. 

Figs.  3.15-3.17  and  3.18-3.20  are  time  sequences  of  the  folding  event;  the  first  se¬ 
quence  shows  the  relationship  between  the  tropopause  fold  and  thermal  field,  and  the 
second  the  relationship  between  the  fold  and  associated  jet.  Perhaps  the  most  important 
feature  of  these  time  series  is  their  general  agreement  with  the  growing  base  of  observa¬ 
tional  data  (e.g.  Shapiro,  1987).  In  these  diagrams  we  see  clearly  the  evolution  from  the 
initial  trough  of  cold  air  shown  in  figs.  3. 3-3. 6,  to  a  lowering  of  the  tropopause,  and  fi¬ 
nally  to  a  narrow,  deepening  front  with  a  statically  stable  layer,  significant  vertical  motion 
along  isentropic  surfaces  through  the  fold,  and  wind  shearing  in  the  vicinity  of  the  front 
and  jet  streak.  At  the  48  hour  point  especially  there  is  a  characteristically  thin  layer  of 
stratospheric  air  within  the  fold  bounded  both  above  and  below  by  tropospheric  air. 

Our  reproduction  of  upper  frontogenesis  processes  in  the  absence  of  an  initial  surface 
temperature  gradient  confirms  previous  results  of,  for  example,  Buzzi  et  al.  (1981).  Fig. 
3.21  shows  perhaps  more  clearly  than  previous  diagrams  the  extent  of  upper  frontogenesis 
obtainable  when  the  deformation  field  varies  vertically.  Among  the  most  notable  features 
of  the  model  is  the  folded  tropopause  associated  with  the  internal  front.  The  realism  of 
the  folding  process,  as  well  as  the  horizontal  wind  shears  developed  near  the  jet  core,  are 
enhanced  by  the  use  of  the  geostrophic  horizontal  coordinate.  Such  increased  resolution 
of  field  parameters  has  been  documented  before  by  Hoskins  and  West  (1979),  for  example, 
but  not  as  applied  in  the  current  results.  Other  recurrent  features  include  the  adiabatic 
exchange  between  troposphere  and  stratosphere  through  the  fold,  and  a  realistic  jet  stream 
in  terms  of  both  proximity  to  the  front  and  intensity.  These  features  taken  together  are 
remarkably  similar  to  observations  of,  for  example,  Shapirc  (1985),  where  an  upper  front, 
containing  primarily  stratospheric  air  is  bounded  both  above  and  below  by  tropospheric 
air.  The  fact  that  nearly  adiabatic  motion  occurs  through  the  model-produced  fold  re¬ 
inforces  the  notion  that  anomalous  lamina  of  ozone  can  be  injected  to  the  troposphere 
through  such  a  tropopause  folding  event  (Danielsen,  1980). 
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Figure  3.15:  Thermal  field  at  24  hours  in  physical  space  with  pressure  (mb)  as  the  vertical 
coordinate.  The  lowest  isotherm  is  285  K  with  increments  of  5  K.  Dashed  line  is  the 
constant  PV  line  2  •  10"7  K/(Pa-s),  representing  the  dynamic  tropopause.  Distance  is  in 
km  on  the  horizontal  axis. 
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Figure  3.16:  Same  as  fig.  3.15,  except  for  36  hours. 
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Figure  3.17:  Same  as  fig.  3.15,  except  for  48  hours. 
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Figure  3.18:  Along-front  geostrophic  winds  at  24  hours  in  physical  space  with  pressure  as 
the  vertical  coordinate.  Dashed  line  is  the  same  as  in  fig.  3.15.  Isotachs,  drawn  every  4 
m/s,  show  the  jet  core  above  the  developing  fold. 
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Figure  3.19:  Same  as  fig.  3.18,  except  for  36  hours. 
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Figure  3.20:  Same  as  fig.  3.18,  except  for  48  hours. 


45 


Chapter  4 


FRONTOGENESIS  BY  GENERALIZED  VERTICAL  WIND  SHEARING 

4.1  Comments 

In  Chapter  3  we  retraced  the  classical  frontogenesis  procedure  of  Hoskins  et  ai.  Using 
a  generalized  geostrophic  deformation  field  which  varied  vertically,  as  well  as  an  improved 
initialization  of  the  potential  vorticity  field,  we  were  able  to  create  a  realistic  tropopause 
folding  event  which  bears  resemblance  to  the  aforementioned  observed  frontal  parameters. 
An  essential  feature  in  the  success  of  our  model  is  the  effective  shearing  nature  of  the 
assumed  wind  field.  Such  an  assumption  is  not  unreasonable  since  significant  vertical  wind 
shearing  is  known  to  exist  in  baroclinic  waves.  Thus  we  can  think  of  our  deformation  field 
as  simply  one  example  of  a  typical  shearing  environment  where  fronts  form,  and  generalize 
our  results  by  eliminating  the  deformation  constraint.  In  its  place  the  wind  will  now  be 
specified  as  a  function  of  the  vertical  coordinate. 

Our  initialization  of  the  a *  field  in  Chapter  3  allowed  for  the  formation  of  surface 
fronts  even  though  the  lower  boundary  was  defined  as  an  isentropic  surface.  Here  we 
examine  another  prescription  of  the  initial  a*  field,  with  emphasis  on  the  elimination  of 
surface  fronts  so  as  to  isolate  internal  processes. 

4.2  Isolating  Internal  Processes 

As  an  example  of  a  more  general  zonal  flow  we  set  ug  =  ug(8),  and  vg  =  v'g  ( x,9,t ) 
where  ug  is  any  function  of  9  and  is  as  defined  in  (3.2).  As  before  we  assume  a 
straight  front  oriented  north-south  with  d/dy  =  0,  and  using  the  definitions  (2.10)  and 
our  generalized  geostrophic  wind  field,  we  have  d/dY  =  0  also.  Equation  (2.8)  reduces  to 
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which  has  as  its  solution 


cr*(X,Q,T)  =  of  (X  —  iisT,0), 


with  <Tj  now  given  by 


tanh 

[■ 

^0x  —  0) 

—  tanh 

a  (0-0)] 

tanh 

a ^0x  —  0) 

-  tanh 

a(0B-0)] 

(<7B  -  ^t) 


+Acr  (  - 4  (9  )  tanh  (/3X) , 


"  V  (0T-0B)2  J  ^ 

and  with  0b  and  0x  as  given  in  Chapter  3,  a b  =  1700  Pa/K,  <rx  =  100  Pa/K,  Act 
half  the  horizontal  variation  of  potential  density  in  the  domain,  a  =  1/50,  0  =  1/<5X, 
n  >  1,  and  0  again  a  constant  global  average  potential  temperature  of  the  tropopause. 
Notice  that  at  the  upper  and  lower  boundaries  (0x  and  0b),  equation  (4.3)  reduces  to 
the  constant  surfaces  <?x  and  <tb-  la  addition,  the  function  (4.3)  has  been  chosen  so  that, 
when  integrated,  the  upper  and  lower  boundaries  are  also  the  constant  pressure  surfaces 
px  and  pa  respectively,  which  satisfies  our  original  boundary  assumptions  in  equations 
(2.26b)  and  (2.26c).  The  associated  diagnostic  equations  for  the  interior  and  boundaries 
are  unchanged  from  Chapter  3: 


a2M* 

OQdX 


+  /2iv*  =  o, 


=  nT  at  0  =  ©x, 


5M*  1  / dM*\ 2 

9¥-Mtv(ir)  =°  at  e  =  0B- 


M*(0,T)  =  M*(0,O)  at  X=-XL,XL. 


(4.4a) 

(4.4b) 

(4.4c) 

(4.4d) 


As  can  be  seen  from  (4.4d),  in  these  experiments  it  was  possible  to  narrow  the  domain 
width  somewhat  to  a  more  manageable  size  (here  we  typically  have  Xl  =  3A).  Once 
the  zonal  geostrophic  winds  are  known  or  prescribed,  our  computational  scheme  is  also 
unchanged:  compute  a"  at  time  T;  use  this  field  to  evaluate  M*;  use  these  values  in 
equation  (2.22)  to  produce  the  wind  and  mass  fields. 
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4.3  Results  and  Comparative  Discussion 

There  are  noteworthy  differences  between  the  experiments  here  and  those  in  Chapter 
3.  By  excluding  the  formation  of  surface  fronts  we  isolate  internal  processes,  but  simul¬ 
taneously  remove  possible  interactions  between  the  two  through  vertical  motion.  Several 
functions  were  considered  as  defining  the  ug  field.  Examples  of  model  runs  which  met  with 
somewhat  less  success  included  a  westerly  jet  core  in  the  upper  troposphere  with  wind 
speed  decreasing  both  above  (where  they  become  easterly  in  the  stratosphere)  and  below 
in  a  symmetric  Gaussian  fashion;  specifically,  retaining  our  notation  from  Chapter  3  we 
set 


In  these  model  runs  no  upper  front  was  produced,  but  the  tropopause  did  fold  in  an 
unrealistic  fashion. 

Perhaps  a  more  realistic  vertical  wind  distribution  is  given  if  the  easterlies  above  the 
front  are  significantly  weaker  than  the  westerlies  in  the  vicinity  of  the  jet  core: 


e-se 


with  90  =  69/ 3,  ug  a  specified  constant,  and  n  any  positive  odd  integer.  Higher  values  of 
n  result  in  a  larger  difference  between  the  maximum  wind  speeds  in  the  troposphere  and 
stratosphere,  respectively. 

We  can  imagine  an  initial  state  where  there  is  a  weak  horizontal  gradient  of  o~  cor¬ 
responding  to  a  low-amplitude  upper  trough.  After  a  period  of  time  during  which  the 
geostrophic  wind  field  (equation  4.5)  is  allowed  to  act,  a  front  will  form  due  to  differential 
advection  by  ug  with  respect  to  0-surfaces.  Results  using  (4.5)  in  (4.1)  are  shown  in  figs. 
4. 1-4.6:  4. 1-4.3  are  a  time  sequence  of  the  potential  density  field,  and  4. 4-4. 6  are  the 
same  sequence  for  the  thermal  field  with  a  representative  tropopause  value  of  potential 
vorticity. 

As  can  be  seen  from  these  figures,  the  upper  front  as  given  by  the  a*  field  is  in  lesser 
agreement  with  observations  than  results  in  Chapter  3  since  each  individual  isoline  has 
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Figure  4.1:  Potential  density  field  in  physical  (x,p)  space  after  ug,  as  given  in  equation 
(4.5),  has  been  allowed  to  act  for  a  period  of  time.  Isolines  are  every  200  Pa/K,  ranging 
from  200  to  1600.  Distance  is  in  km  on  the  horizontal  axis. 
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Figure  4.2:  Same  as  fig.  4.1,  except  at  a  later  time. 


Figure  4.3:  Same  as  fig.  4.2,  except  at  a  later  time. 
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Figure  4.4:  0-field  (K)  and  tropopause  {q  =  7  •  10-8  K/(Pa-s))  for  the  same  time  as  fig. 
4.1.  Coordinate  system  is  as  in  Fig.  4.1,  and  isentropes  are  every  5  K  beginning  with  285 
K. 
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Figure  4.5:  Same  as  fig.  4.4  but  for  the  same  time  as  fig.  4.2. 


-2000 


0 


2000 


4000 


Figure  4.6:  Same  as  fig.  4.4  but  for  the  same  time  as  fig.  4.3. 
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a  relatively  shallow  vertical  extent;  that  is,  there  is  no  single  value  of  q  defining  a  deep 
upper-front.  Taken  as  a  whole,  however,  the  folding  process  extends  deeper  into  the  tro¬ 
posphere  than  our  earlier  results.  The  thermal  field  here  is  again  in  rough  agreement  with 
observations,  especially  with  respect  to  its  relationship  with  the  PV  field.  Although  the 
meridional  wind  contained  realistic  magnitude  and  shears,  its  proximity  to  the  front  was 
not  well  correlated.  Also,  the  initial  winds  were  quite  strong  and  the  jet  did  not  develop 
with  the  developing  front — these  factors  are  attributable  to  the  implicit  relationship  be¬ 
tween  our  assumed  erf  field  and  the  associated  M*  field.  As  stated  earlier  equation  (4.3) 
has  the  advantage  of  producing  a  fiat  PV  field  on  the  upper  and  lower  boundaries.  A 
drawback,  however,  is  that  this  particular  function  would  not  allow  a  sharp  vertical  gra¬ 
dient  near  the  tropopause  while  still  allowing  for  realism  in  other  respects,  such  a  defining 
a  monotonically  decreasing  function  of  0. 


Chapter  5 


SUMMARY  AND  CONCLUSIONS 

We  have  sought  to  bridge  a  gap  in  understanding  the  formation  of  internal  fronts 
within  a  framework  which  considers  them  as  a  byproduct  of  baroclinic  waves.  Building 
on  the  accomplishments  of  Hoskins  et  al.  by  using  geostrophic/isentropic  coordinates 
with  a  geostrophic  deformation  field  as  a  forcing  mechanism,  the  objective  of  simulating 
tropopause  folding  in  a  balanced  model  much  as  Keyser  ef  al.  have  done  in  a  primitive 
equation  model  has  largely  been  successful.  We  also  departed  from  Keyser’s  work  in  that 
horizontal  shear  is  neglected  for  our  two-dimensional  case;  instead  we  considered  the  effects 
of  vertical  wind  shears  growing  out  of  the  baroclinic  wave  accompanied  by  a  continuous 
potential  vorticity  field.  The  choice  of  coordinates  simplifies  the  dynamics  by  implicitly 
including  both  ageostrophic  and  vertical  motions  in  a  system  of  equations  which  are  as 
simple  as  the  quasigeostrophic  equations,  yet  are  applicable  in  a  broad  and  important  class 
of  motions  (high  relative  vorticity  but  not  highly  curved  flow)  where  quasi-geostrophic 
theory  breaks  down. 

Taken  as  a  whole  our  results  show  that  internal  fronts  can  form  in  a  balanced  two- 
dimensional  model  for  two  distinct  cases — in  the  absence  of  a  surface  temperature  gradient, 
where  both  surface  and  upper  fronts  can  form,  and  when  surface  fronts  are  eliminated 
as  a  possibility — when  the  working  definition  of  the  tropopause  is  a  surface  of  constant 
potential  vorticity.  The  first  case,  which  is  covered  in  Chapter  3,  provides  the  best  results. 
With  a  horizontal  temperature  gradient  of  25  K  across  the  upper  trough,  the  model 
produced  a  fold  extending  down  to  approximately  the  500  mb  level.  A  jet  of  nearly 
100  kt  developed  along  the  front  and  above  the  fold  with  realistic  wind  shears  when  the 
domain  was  shifted  from  geostrophic  back  to  physical  space.  Isentropic  surfaces  pass 
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through  the  front,  showing  that  nearly  adiabatic  motions  of  dry  stratospheric  air  can 
be  advected  into  the  mid  troposphere.  Thus,  there  is  a  correlation  between  the  thermal, 
wind  and  mass  fields  which  strongly  resembles  previously  observed  features  of  upper  fronts. 
These  observations  used  various  combinations  of  aircraft,  dropwindsonde,  and  rawindsonde 
data.  The  second  case  (Chapter  4)  reinforces  some  results  of  Chapter  3:  a  fold/upper 
front  results  from  the  model,  and  its  relationship  to  the  potential  temperature  field  shows 
similarities  as  well.  Parameters  near  the  front  were  more  highly  resolved,  also,  again  due 
to  the  coordinate  shift. 

Based  on  our  results  there  seems  to  be  sufficient  evidence  to  continue  the  study 
frontogenesis  using  a  simple,  two-dimensional,  balanced  model  approach.  Here  we  were 
able  to  simplify  the  dynamics  of  the  system  to  essentially  two  equations — one  predictive 
and  one  diagnostic.  Due  to  the  geostrophic  momentum  approximation,  these  equations 
have  a  more  general  applicability  than  quasigeostrophic  theory.  By  allowing  vertical  wind 
shear  to  be  part  of  the  initialization  process,  a  trough  of  cold  air  aloft  bordered  by  warmer 
air  on  either  side  formed  a  rather  deep  upper  front  which  was  a  result  of  a  tropopause 
folding  event.  Unlike  previous  results  using  balanced  models,  reproduction  of  the  true  fold 
was  not  dependent  upon  the  shift  back  to  physical  space — use  of  geostrophic  coordinates 
simply  enhanced  the  realism  of  the  final  product. 

The  introduction  of  more  generalizations  into  the  model  leaves  more  areas  which 
would  benefit  from  further  study.  First  it  would  be  useful  to  allow  isentropic  surfaces  to 
intersect  the  lower  boundary,  possibly  by  utilizing  the  positive  definite  scheme  of  Arakawa. 
Results  would  be  more  general  and  realistic  of  observed  surface  features  for  most  cases 
of  upper  frontogenesis.  Accomplishing  this  alone  would  be  noteworthy  by  virtue  of  its 
application  to  other  areas  since  the  major  drawback  of  how  to  handle  the  lower  boundary 
condition  would  be  solved.  At  the  same  time,  the  benefits  of  the  isentropic  vertical  co¬ 
ordinate  which  have  made  it  an  increasingly  popular  choice,  such  as  retention  of  vertical 
motion  for  the  adiabatic  case,  would  be  retained.  Second,  more  can  be  done  to  correlate 
observational  data  with  our  results.  For  example,  vertical  temperature  profiles  through 
a  model  produced  fold  may  be  similar  to  observations  of  Danielson  (1980).  Introducing 
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moisture  and  ozone  as  advected  quantities  should  reveal  a  dry  tropospheric  layer  with 
both  anomalously  high  ozone  values  and  low  liquid  water  content  at  the  level  of  the  fold, 
indicative  of  the  aforementioned  exchange  process  between  stratosphere  and  troposphere. 
Again,  this  would  correlate  well  with  observations.  Third,  although  we  included  the 
method  for  calculating  the  secondary  circulation  across  the  front,  recovering  these  fields 
from  the  model  would  provide  insight  to  the  tilting  and  subsidence  effects  they  produce. 
Finally,  more  work  can  be  done  defining  the  initial  fields  for  the  Chapter  4  case  (where  the 
upper/lower  boundaries  are  constant  potential  vorticity  surfaces)  so  that  there  is  initially 
a  weaker  geostrophic  flow  along  the  front,  and  the  resulting  jet  which  develops  over  time 
is  better  correlated  with  the  formation  of  the  front. 

In  retrospect  it  seems  likely  that  the  major  results  contained  herein  could  also  be 
achieved  for  a  two-layer  model.  Such  an  application  would,  by  virtue  of  its  link  to  both 
previous  work  and  current  results,  help  in  advancing  the  idea  that  much  remains  to  be 
learned  about  the  subject  by  employing  the  two-dimensional  balanced  approach. 
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